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ABSTRACT 

Precise radial velocity measurements from HIRES on the Keck I telescope are pre- 
sented for 40 stars in the outer halo globular cluster NGC 2419. These data are used 
to probe the cluster's stellar mass function and search for the presence of dark mat- 
ter in this cluster. NGC 2419 is one of the best Galactic globular clusters for such a 
study due to its long relaxation time (T r o ~ 10 10 yr) and large Galactocentric dis- 
tance (Rgc ~ 90 kpc) — properties that make significant evolutionary changes in the 
low-mass end of the cluster mass function unlikely. We find a mean cluster velocity 
of <v r >= —20.3 ± 0.7 km/sec and an internal velocity dispersion of a = AAA ± 0.48 
km/sec, leading to a total mass of (9.0±2.2) -10 5 M and a global mass-to-light ratio of 
MjLy — 2.05 ±0.50 in solar units. This mass-to-light ratio is in good agreement with 
what one would expect for a pure stellar system following a standard mass function at 
the mctallicity of NGC 2419. In addition, the mass-to-light ratio does not appear to 
rise towards the outer parts of the cluster. Our measurements therefore rule out the 
presence of a dark matter halo with mass larger than ~ 10 7 M Q inside the central 500 
pc, which is lower than what is found for the central dark matter densities of dSph 
galaxies. We also discuss the relevance of our measurements for alternative gravita- 
tional theories such as MOND, and for possible formation scenarios of ultra-compact 
dwarf galaxies. 
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1 INTRODUCTION 

An important unanswered question in modern astrophysics 
concerns the mechanism(s) by which stars and star clusters 
form. This is especially true in the case of globular clusters, 
which, due to their large ages, are probes of the evolution 
and assembly history of galaxies in the early universe. 

Various theories have been proposed concerning the for- 
mation process of globular clusters. A number of authors 
proposed that globular clusters formed at the centers of 
dark matter halos (e.g. Peebles 1984, West 1993, Cen 2001, 
Kravtsov & Gnedin 2005). Such a cosmologically-motivated 
formation scenario for globular clusters is indicated by the 
extreme age of most clusters, their low metallicities, and 
their characteristic masses (which are close to the Jeans 
mass at recombination). The fact that most clusters in the 
Milky Way contain no measurable amount of dark matter is 
not in contradiction to these theories, since, over the course 
of a Hubble time, dark matter would get depleted from the 
central regions of globular clusters due to the dynamical 



friction of stars against the much lighter dark matter parti- 
cles (Baumgardt & Mieske 2008). As a result, dark matter 
tends to get pushed towards the outer cluster parts, where 
it is much more difficult to detect and can be more eas- 
ily stripped away by the tidal field of the Milky Way (e.g., 
Mashchenko & Sills 2005). 

A significant dark matter content might also be the ex- 
planation for the high mass-to-light ratios of ultra-compact 
dwarf galaxies (UCDs), which have been discovered in spec- 
troscopic surveys of nearby galaxy clusters since the late 
1990s (Hilker et al. 1999; Drinkwater et al. 2000). UCDs arc 
bright (— 11 < Mv < —13.5 mag) and compact (7 < < 
100 pc) stellar systems that: (1) appear to obey a set of 
structural scaling relations distinct from those of globular 
clusters; and (2) have mass-to-light ratios that are, on aver- 
age, about twice as large than those of globular clusters of 
comparable metallicity, and somewhat larger than what one 
would expect based on simple stellar evolution models as- 
suming a standard stellar initial mass function (Hasegan et 



2 Baumgardt et al. 



al. 2005; Dabringhausen, Hilker & Kroupa 2008; Mieske et 
al. 2008). One possible formation scenario for UCDs could be 
adiabatic gas infall into the center of a dwarf galaxy, which 
also funnels dark matter into the center (Goerdt et al. 2008). 
Later tidal disruption of these galaxies would leave only the 
central star cluster behind. Since some of the most mas- 
sive Galactic globular clusters have masses and sizes only 
slightly smaller than UCDs (e.g. Rejkuba et al. 2007), it is 
conceivable that some of them may have formed in a similar 
way. 

On the other hand, observations of interacting and star- 
burst galaxies have shown that star clusters with masses 
comparable to globular clusters also form during major 
mergers between galaxies (Whitmore & Schweizer 1995; 
Whitmore et al. 1999; Goudfrooij et al. 2001) in the lo- 
cal universe, which makes it possible that globular clusters 
formed by similar processes in the early universe (e.g. Searle 
& Zinn 1978, Ashman & Zepf 1992). In such a scenario, glob- 
ular cluster formation is driven mainly by gas-dynamical 
processes and the globular clusters should not contain sig- 
nificant amounts of dark matter. 

Globular clusters in the outer halo of the Milky Way 
are ideal objects to test for the presence of dark matter 
and therefore examine the formation mechanism of globular 
clusters. Because the tidal field is weak in the remote halo, 
tidal stripping of the outer cluster regions will be much less 
important than for globular clusters orbiting close to the 
Galactic center. Baumgardt & Mieske (2008) have shown 
that, due to dynamical friction, dark matter gets depleted 
from the central regions of a star cluster on a timescale 
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where Mmt is the total cluster mass, is the half-mass 
radius and m is the average stellar mass in the cluster. For 
clusters with Tfhc < Tkubbie, dark matter will be removed 
from the central parts within a Hubble time. 

Fig. 1 shows the dynamical friction timescales of Galac- 
tic globular clusters. Here the cluster masses and half-mass 
radii, ru, were calculated from the absolute magnitudes, 
projected half-mass radii rn v and distances given by Harris 
(1996), assuming M/L = 2.5, which is appropriate for most 
metal-poor clusters with [Fe/H] < —1.2 and r% = 1.33rft P , 
which is approximately correct for most King models. It can 
be seen that the majority of galactic globular clusters have 
friction timescales much less than a Hubble time. According 
to the results of Baumgardt & Mieske (2008), dark mat- 
ter would exist in such clusters only in the outer parts, 
where it is difficult to detect due to low stellar densities 
and the influence of the Galactic tidal field on the stellar 
velocities. Only a few clusters have friction timescales sig- 
nificantly longer than a Hubble time and are therefore good 
candidates to look for primordial dark matter. Among the 
clusters with long dynamical friction timescales are Pal 3, 
Pal 4 and Pal 14, which are currently being investigated by 
Jordi et al. (2009) in an effort to test alternative gravity 
theories like MOND (Milgrom 1983a,b). 

The Galactic globular cluster with the longest dynami- 
cal friction time, and therefore the smallest amount of dark 
matter depletion, is NGC 2419. With a Galactocentric dis- 
tance of Rgc ~ 90 kpc, NGC 2419 is also one of the most 
remote Galactic globular clusters, meaning that tidal strip- 
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Figure 1. Dynamical friction time, Tp r i c , of stars against lighter 
dark matter particles for Galactic globular clusters. Most globu- 
lar clusters have friction timescales of less than a Hubble time, 
meaning that dark matter would have been depleted from their 
centers if they formed as a mix of dark matter and stars. Only 
a few extended clusters have friction times longer than a Hubble 
time and should therefore still retain dark matter in their cen- 
ters. With the longest friction timescale of all Galactic globular 
clusters, NGC 2419 is a promising target for a search for dark 
matter. 



ping was likely to be least effective in this cluster. Finally, 
NGC 2419, as one of the most massive Galactic globular 
clusters (see Fig.l), represents a possible Local Group ana- 
logue to a UCD. Its half-light radius of r hp = 47.5" « 19 pc 1 
(see §3.1) is also ~ 6x larger than is typical for globular clus- 
ters in the Milky Way and external galaxies (Jordan et al. 
2005), although still a factor of ~ 2 smaller than the most 
extreme UCDs (see, e.g., Table 5 of Mieske et al. 2008). 
For all of these reasons, NGC 2419 is an excellent place to 
search for the presence of dark matter and thereby test dif- 
ferent formation scenarios of globular clusters and, possibly, 
UCDs. 

The only published velocity dispersion for NGC 2419 is 
that of Olszewski, Pryor & Schommer (1993) who, based on 
MMT radial velocity measurements with a median precision 
of 1.6 km/sec for 12 stars, found a mean radial velocity of 
v c = —20.0 ± 0.9 km/sec and an intrinsic velocity disper- 
sion of cr c = 2.7 ± 0.8 km/sec. Fitting single-mass, isotropic 



1 At our adopted distance for NGC 2419, (m- M) = 19.60 mag 
(Ripepi et al. 2007), l"corresponds to 0.40 pc. 
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King (1962) models to their radial velocities and the clus- 
ter surface brightness profile yielded a mass-to-light ratio of 
M/L v ~ 0.7 ± 0.4 Mq/L v ,q — one of the smallest values 
measured for a Galactic globular cluster and much lower 
than the M/Lvs found for some UCDs. 

In this paper, we report on radial velocity measurements 
for 44 candidate red giant branch (RGB) stars in the direc- 
tion of NGC 2419 taken with the High Resolution Echelle 
Spectrometer (HIRES; Vogt et al. 1994) on the Keck I tele- 
scope. The paper is organized as follows: In § 2 we describe 
the details of the observations and data reduction. In § 3 we 
determine the total cluster luminosity, velocity dispersion 
and mass-to-light ratio of NGC 2419, and in § 4 we draw 
our conclusions. 



2 OBSERVATIONS AND DATA REDUCTION 

Candidate RGB stars in NGC 2419 were selected for obser- 
vation with HIRES from the photographic study of Racine 
& Harris (1975), as well as from unpublished CCD photom- 
etry of Peter Stetson. Fig. 2 presents a finding chart for the 
candidate RGB stars that were observed for this program. 

HIRES spectra were obtained on four separate nights 
during two observing runs on the Keck I telescope. The dates 
of the two runs were February [10-11] and March [9-11] 1999 
(i.e., five nights in total). During each run, we limited the 
entrance aperture using the CI decker (0.866 x 7.0") and 
binned the detector 1x2 (i.e., in the spatial direction) to 
reduce the read noise. The spectral resolution for this in- 
strumental configuration is 1Z = 45000. A single readout 
amplifier was used (gain = 2.4 ADU~ ) with the red colli- 
mator and a cross-disperser in first order. Thorium-Argon 
comparison lamp spectra were acquired frequently during 
each night. In all, spectra were acquired for 44 different RGB 
candidates during the two runs, with exposure times rang- 
ing from 180 to 600 sec. Approximately two dozen high-S/N 
spectra for seven distinct IAU radial velocity standard stars 
were also obtained during the course of the program. 

Spectra for program objects and standard stars were 
reduced in an identical manner following the general proce- 
dures described in Cote et al. (1999, 2002). Briefly, the ra- 
dial velocity of each RGB candidate was measured by cross- 
correlating its spectrum against that of a master template 
created during each run from the observations of IAU stan- 
dard stars. In order to minimize possible systematic effects, 
a master template for each observing run was derived from 
an identical subsample of IAU standard stars. From each 
cross-correlation function, we measured both v r , the helio- 
centric radial velocity, and Rtd, the Tonry & Davis (1979) 
estimator of the strength of the cross-correlation peak. 

The NGC 2419 observations discussed here were taken 
as part of a broader program to study the internal dynam- 
ics of outer halo globular clusters, some results of which 
have been reported elsewhere (i.e., Cote et al. 2002, Jordi et 
al. 2009). During this program, which spanned seven Keck 
observing runs in 1998 and 1999 (13 nights in total), we ob- 
tained 53 distinct radial velocity measurements for 23 differ- 
ent RGB and subgiant stars belonging to Palomar 3, Palo- 
mar 4, Palomar 5, Palomar 14, and NGC 7492, in addition 
to NGC 2419. Using these repeat measurements and fol- 
lowing the procedures described in Vogt et al. (1995), we 



derived an empirical relationship between our radial veloc- 
ity uncertainties, e(v r ), and the strength of cross correlation 
peak: e(v r ) = a/(l + R TD ), where a = 9.0±ll km/sec (90% 
confidence limits). 

Our measurements are summarized in Table 1. The 
name of each program star is recorded in Column (1), with 
most identifications taken directly from Racine & Harris 
(1975). The exception is for stars with IDs beginning with an 
' S ' prefix, which denotes targets selected from our (unpub- 
lished) photometric catalogue. The second column of this ta- 
ble gives the ID from Olszewski, Pryor & Schommer (1993) if 
the star appeared in the earlier MMT study. Columns (3)- 
(9) record the position in right ascension and declination, 
the radial distance from the cluster center, V magnitude and 
(B-V) colour from our catalogue, the number of independent 
HIRES radial velocity measurements, and the adopted he- 
liocentric radial velocity, v r . For stars with multiple velocity 
measurements, the velocity reported in Column (9) refers 
to the weighted mean of the individual measurements. We 
assumed that the cluster center is right ascension (J2000) 07 
38 08.47 and declination +38 52 56.8 in order to calculate 
radial distances of stars. 

In the left panel of Fig. 3 we show a V,(B-V) 
colour-magnitude diagram (CMD) for stars within 250"of 
NGC 2419 — extending from the tip of the RGB to the lower 
subgiant branch — that shows the location of our program 
stars in the CMD. For comparison, the solid curve shows an 
isochrone from Dotter et al. (2008) having T = 12.3 Gyr, 
and matching the chemical abundances of the cluster: e.g., 
[Fe/H] = -2.32 dex and [a/Fe] = +0.2 dex (e.g., Shetrone, 
Cote & Sargent 2001). In making this comparison, we have 
adopted a reddening of E(B — V) = 0.10 mag — intermedi- 
ate between the values of 0.11 mag and 0.08 mag reported 
by Harris et al. (1997) and Ripepi et al. (2007) — and a 
de-reddened distance modulus of (m — M) = 19.60 mag 
(Ripepi et al. 2007). Note that all of our program stars, and 
those from the previous MMT program, are located within 
~ 2 mag of the RGB tip. The right panel of Fig. 3 shows a 
magnified view of this part of the CMD. The dashed lines 
show a region of ±0.15 mag width centered on the isochrone 
— an interval chosen to roughly match the observed width of 
the RBG sequence and encompass all of the probable radial 
velocity members. 

Since NGC 2419 is located at a Galactic latitude of 
b ~ 25°, it is important to consider the possibility of con- 
tamination by foreground disk stars. To gauge the extent of 
such contamination, we show in the right panel of Fig. 3 the 
results of one simulation using the Besangon Galaxy model 
(Robin et al. 2003), for field stars in the range 0-100 kpc 
and located within a 500" x 500" region in the direction of 
NGC 2419. This simulation yields a total of 10 stars hav- 
ing magnitudes and colours that would place them within 
the region bounded by the dashed lines in the right panel 
of Fig. 3. However, as we shall see in §3.2, only two stars 
in this simulation have radial velocities within ±3a of the 
cluster mean. Since our radial velocity survey is by no means 
complete within this simulated region, we conclude that our 
final sample should have minimal foreground contamination 
with 1 interloper, at most, expected within ±3<r of the clus- 
ter's systemic velocity. The final column of Table 1 gives our 
division of the sample into cluster members and foreground 
stars. In most cases, the assignments are unambiguous, al- 
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Table 1. Radial velocities for candidate red giants in NGC 2419 



ID 


093 




R.A. 


Decl. 


R 


V 


(B-V) 


N Vr 


v r 


Member? 






(J2000.0) 


(J2000.0) 


(arcsec) 


(mag) 


(mag) 




(km s" 1 ) 




(1) 


(2) 




(3) 




(4) 




(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


III89 


III89 


07 


38 


09.79 


+38 


50 


45.4 


132.3 


17.25 


1.42 




Of. n^-r-fl Q9 


Y 


Do 




07 


38 


01.92 


+38 


53 


15.4 


78 7 


1 7 9fi 


1 44 




ofi ^9-1-0 fil 


Y 


S2 




07 


38 


09.90 


+38 


52 


00.7 


58.5 


17.27 


1.55 


I 


- J. U.UUIU.OO 


Y 


S4 




07 


38 


07.71 


+38 


53 


15.7 


20.9 


17.27 


1.51 


I 


oa 97-i-fl fifi 


Y 


S5 




07 


38 


10.07 


+38 


52 


54.3 


18.9 


17.29 


1.48 


I 


i f. 41 _|_n fj 


Y 


S12 




07 


38 


06.00 


+38 


53 


37.5 


49.9 


17.40 


1.39 


I 


99 1 ftJ-fl f.Q 


Y 


S6 




07 


38 


06.89 


+38 


53 


33.5 


41.1 


17.41 


1.28 


I 


1 f. 9A-r-n 71 


Y 


S14 




07 


38 


11.58 


+38 


53 


05.8 


37.4 


17.42 


1.32 


I 


1 9 99-Un fil 


Y 


S9 




07 


38 


10.18 


+38 


53 


09.6 


23.7 


17.42 


1.24 


I 




Y 


*y± 




07 


38 


04.90 


+38 


51 


50.6 


78 9 


1 / .4:0 


l.OJ 


^ 


90 Qfi-l-1 98 


Y 


S10 




07 


38 


16.92 


+38 


53 


35.0 


105.8 


17.43 


1.35 




90 7^-1-0 ^8 


Y 


S17 




07 


38 


08.01 


+38 


52 


53.1 


6.5 


17.46 


1.22 




1 1 9-r-O 81 
_LU . J. ZZC U.Ol 


Y: 


S15 




07 


38 


09.69 


+38 


52 


40.9 


21.4 


17.48 


1.43 


I 


1 7Q-I-0 ^8 
-_LU. 1 yiu.uo 


Y 


S20 




07 


38 


07.60 


+38 


53 


34.1 


38.7 


17.50 


1.38 


I 


00 fifi-i-o ^ 


Y 


S16 




07 


38 


14.21 


+38 


52 


36.0 


70.2 


17.51 


1.36 


I 


_1 Q QQ+n Af. 


Y 


S22 


C37 


07 


38 


08.80 


+38 


52 


42.9 


14.4 


17.51 


1.31 


I 


-00 . yoiu .urt 


Y: 


S23 


C80 


07 


38 


09.57 


+38 


53 


38.7 


43.8 


17.54 


1.32 


I 


-99 81 +0 "SS 


Y 


S26 




07 


38 


05.90 


+38 


52 


52.0 


30.4 


17.57 


1.38 


I 


97 1 7-1-0 ^3 


Y 


10 


10 


07 


38 


16.92 


+38 


53 


35.0 


105.8 


17.61 


1.28 


2 




Y 


38 




07 


38 


02.87 


+38 


52 


24.1 


73.1 


17.61 


1.41 


3 


-ZO.UOIU.1U 


Y 


3 




07 


38 


03.48 


+38 


54 


12.1 


95.2 


17.63 


1.34 


I 


90 7(S-|-0 7Q 


N 


S31 




07 


38 


08.64 


+38 


53 


00.0 


3.8 


17.64 


1.19 


I 


9Q RQ-l-0 43 
-ZO .Ut>I U.IO 


Y 


S38 


38 


07 


38 


07.67 


+38 


52 


45.5 


14.7 


17.75 


1.18 


i 


-1 99+n 43 


Y: 


S41 




07 


38 


06.32 


+38 


52 


38.7 


31.0 


17.79 


1.27 


I 


90 97-I-O 44 


Y 


IV93 




07 


37 


56.45 


+38 


51 


22.8 


169.0 


17.80 


1.29 


I 


-90 (S4+f) 6Q 


Y 


1123 




07 


38 


09.58 


+38 


56 


30.0 


213.6 


17.95 


1.14 


I 


90 (S7-|-n 8Q 
zu.u 1 xu.oy 


Y 


4 




07 


38 


05.94 


+38 


53 


50.8 


61.6 


17.97 


1.27 


2 


90 4Q-I-0 83 

-ZiO.iyiU.OO 


Y 


31 




07 


38 


11.15 


+38 


51 


53.4 


70.7 


18.07 


1.17 




1 q 09-1-1 fl8 


Y 


193 




07 


38 


01.00 


+38 


54 


57.8 


149.1 


18.08 


1.54 


I 


70 fil -r-0 QQ 
1 u.unu.yy 


N 


II118 




07 


38 


16.93 


+38 


54 


25.2 


132.5 


18.09 


1.15 


I 


90 7^-1-1 D4 

ZU. 1 UDI J. .U1 


Y 


148 




07 


37 


58.06 


+38 


55 


22.2 


189.5 


18.13 


1.12 


2 


1 q 47-1-0 71 


Y 


1105 




07 


37 


54.13 


+38 


54 


32.7 


192.9 


18.16 


1.16 


1 


-21.43+0.97 


Y 


III86 




07 


38 


09.54 


+38 


50 


31.7 


145.6 


18.23 


1.19 


1 


-18.89+0.99 


Y 


11 




07 


38 


17.81 


+38 


54 


06.1 


129.2 


18.23 


1.10 


1 


48.83+0.64 


N 


III37 




07 


38 


13.72 


+38 


50 


00.5 


186.7 


18.29 


1.05 


1 


-21.30+1.20 


Y 


36 




07 


38 


03.29 


+38 


51 


57.2 


84.9 


18.31 


1.23 


1 


-20.90+1.29 


Y 


35 




07 


38 


03.55 


+38 


51 


51.9 


86.7 


18.32 


1.18 


1 


-18.45+0.98 


Y 


III107 




07 


38 


15.55 


+38 


51 


13.5 


132.3 


18.39 


1.06 


1 


-17.43+0.99 


Y 


1143 




07 


37 


59.55 


+38 


54 


18.2 


132.2 


18.39 


1.39 


1 


45.81+0.62 


N 


12 




07 


38 


18.63 


+38 


53 


52.0 


130.8 


18.60 


1.02 


1 


-17.70+1.16 


Y 


1 




07 


38 


02.83 


+38 


54 


01.7 


92.5 


18.61 


1.01 


1 


-21.63+1.31 


Y 


24 




07 


38 


18.18 


+38 


52 


15.9 


120.5 


18.70 


0.98 


1 


-21.68+1.11 


Y 


1159 




07 


38 


13.97 


+38 


55 


26.7 


163.1 


18.80 


0.98 


1 


-23.20+1.11 


Y 


7 




07 


38 


11.93 


+38 


54 


19.4 


92.0 


19.00 


0.95 


1 


-22.40+1.26 


Y 



though three stars (S17, S22 and S38) are more problemat- 
ical. We shall return to this issue in §3. 



3 RESULTS 

3.1 Surface density profile 

The surface brightness profile of NGC 2419 has most re- 
cently been determined by Bellazzini (2007) using surface 
photometry and star counts from HST/ACS and SDSS. In 
this paper we use his data for our analysis. Fig. 4 shows 
the V-band surface brightness profile of NGC 2419 as deter- 
mined by Bellazzini (2007) and compares it with several dif- 



ferent model profiles. Bellazzini (2007) found that the best- 
fitting King (1962) profile has a central core radius of r c — 
19.2" and a concentration index of c = log(rt/r c ) = 1.35. 
We show this profile by the blue dashed curve in Fig. 4. As 
can be seen, a King profile provides a satisfactory fit to the 
observed surface density profile out to « 300". Beyond this 
point, however, the observed brightness profile declines more 
slowly than the fitted King profile. Such an excess could be 
the signature of extra-tidal stars (e.g., Leon et al. 2000). 
However, at the distance of NGC2419, 300" corresponds to a 
physical distance of about 120 pc. The expected tidal radius 
of NGC 2419 at its current distance is r t ~ 750 pc: i.e., more 
than 5 times as large, if we assume a circular cluster orbit 
and a logarithmic potential with circular velocity v c = 200 
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Figure 2. V-band finding chart for NGC 2419 showing the loca- 
tion of candidate RGB stars with HIRES radial velocity measure- 
ments listed in Table 1. Open circles and squares indicate clus- 
ter members and non-members, respectively. Arrows show the 12 
stars with MMT radial velocity measurements from Olszewski, 
Pryor & Schommer (1993). The cross shows the adopted cluster 
center. North is up and east is to the left in this image, which 
measures 7.'6 X 7.'6. 

km/sec for the Milky Way and a cluster mass of 10 6 Mq. In 
addition, tidal radii often correspond to breaks in the surface 
density distribution but no such break is visible in Fig. 4. 
Hence, unless the cluster orbit is highly eccentric (e > 0.9), 
the outermost stars are most likely bound to NGC 2419 so 
that the tidal radius has not yet been reached by the ob- 
servations, something which was also noted by Ripepi et al. 
(2007). 

We have therefore also tried to fit the surface bright- 
ness profile with Sersic (1968) profiles, which are known to 
provide accurate representations of the surface brightness 
profiles of most early-type galaxies, including UCDs, over 
nearly their entire profiles (e.g., Graham & Guzman 2003, 
Ferrarese et al. 2006, Cote et al. 2007, Evstigneeva et al. 
2007, Hilker et al. 2007). We find that Sersic models that 
best fit the data over the whole range (red dotted curve) 
and those fitted to the data for r > 20" both give good fits 
in the outer profiles but significantly overpredict the central 
surface brightness. In order to obtain a model that fits the 
surface brightness profile over the whole radial range, we 
add a constant density core to the Sersic profile fitted to the 
data for r > 20": 

E(r) = E exp(-(1.9992n - 0.3271) • {r'/r e ) 1/n - 1.0) (2) 



Figure 3. Left panel: Colour-magnitude diagram for NGC 2419 
from our own unpublished photometry, showing the location of 
our program stars. Filled blue circles show the 40 probable mem- 
bers according to our HIRES radial velocities. Open circles show 
the 12 stars with radial velocities measurements from Olszewski 
et al. (1993); there are five stars in common between the two stud- 
ies. Red triangles show the four obvious non-members stars from 
Table 1. The smooth curve shows an 12.3 Gyr isochrone (Dotter 
et al. 2008) with [Fe/H] = -2.32 and [a/Fe] = +0.2 dex, shifted 
by E(B -V) = 0.10 mag and (m - M) v = (m - M) + A v = 
19.91 mag to match the RGB. Right panel: Magnified view of the 
dotted region from the previous panel. The dashed curves show 
a ±0.15 mag region centered on the isochrone. The red asterisks 
illustrate one simulation of the expected contamination in region 
of the colour-magnitude diagram from Galactic foreground stars 
in the direction of NGC 2419 using the Besangon Galaxy model 
(Robin et al. 2003). See text for details. 

where r' = \/r 2 + r§. A \ 2 minimization gives 14" as best 
fitting value for r c . The solid red curve in Fig. 4 shows the 
surface brightness profile of the outer Sersic profile once a 
constant density core with r c — 14" has been added. With 
this model we obtain a very good fit to the overall surface 
brightness profile. The lower panel shows the residuals be- 
tween our fitted model and the observed data. At most radii, 
the differences between the observed surface density and our 
model are within the observational uncertainties. 

Table 2 summarises our results. For a King model fit, we 
find that the best-fitting profile has r c — 19" and c = 1.39, 
which agrees reasonably well with Bellazzini (2007). The 
other rows list the parameters of our best-fitting Sersic pro- 
files. The integrated V magnitudes are calculated separately 
for each surface density profile. They are, in general, some- 
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Figure 4. Upper panel: V-band surface brightness profile of 
NGC 2419 as determined by Bellazzini (2007). The blue, short- 
dashed curve shows the best fitting single-mass, isotropic King 
(1962) model. This model cannot fit the "excess" light at large 
radii. The dotted and long-dashed red curves are fits of two Sersic 
profiles that fit the outer parts of the cluster well but not the cen- 
tral region. The solid curve is a Sersic profile with an added core 
of size r c = 14". It provides a good match to the density profile 
and is within the reported observational uncertainties at most 
radii (bottom panel). 

what smaller than the values by Bellazzini (2007) , who found 
V=10.47 ± 0.07 mag. Because the King profile underpre- 
dicts the light in the outer parts while the two Sersic pro- 
files significantly overpredict the light in the inner parts, 
we consider the values from the cored Sersic profile to be 
most reliable and we will use them in the remainder of 
this paper. For the cored Sersic profile, we obtain a to- 
tal luminosity of V=10.57 mag and a projected half-light 
radius of r^ v = 47.5". The latter value agrees quite well 
with the half-light radius found by Trager, King & Djor- 
govski (1995), rhp = 45.7". With a distance modulus of 
(m — M)o = 19.60 ± 0.05 mag and a reddening of E(B- 
V)=0.08 mag (Ripepi et al. 2007), these values lead to a 
total absolute magnitude and projected half-mass radius of 
NGC 2419 of M v = -9.28 mag and r hp = 19.2 pc respec- 
tively. 

3.2 Velocity dispersion of NGC 2419 

Table 1 lists 44 stars in the direction of NGC 2419 with ra- 
dial velocities measured in the course of this study. Rejecting 
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the four stars with the most discrepant velocities and conse- 
quently identified as non- members in Table 1 (e.g., S3, 193, 
11 and 1143) leaves us with a sample of 40 stars having radial 
velocities in the range —35 ^ v r ^ —5 km/sec. Although this 
sample shows a clear peak at ~ —20 km/sec corresponding 
to the cluster systemic velocity, there are three stars that 
differ by ~ 10 km/sec or more from the apparent cluster 
mean. The upper panel of Fig. 5, which plots radial velocity 
against distance from the cluster center, reveals all three to 
be among the centrally concentrated stars (i.e., with radii of 
r < 15") and thus unlikely to be interlopers. Nevertheless, 
we now pause to consider the extent to which our sample 
could be contaminated by foreground disk stars. 

As noted in §2, one simulation of the expected Galactic 
foreground using the Besangon model (Robin et al. 2003) 
suggests we might expect a total of ~ 10 interloping field 
stars: (1) within a 500" x 500" region along this line of sight; 
and (2) falling along the cluster RGB shown in the right 
panel of Fig. 3. In the lower panel of Fig. 5 we compare 
the observed velocity distribution of radial velocities (open 
blue histogram) with that from the simulation (filled red 
histogram). The foreground stars plotted here represent all 
stars in this velocity range with magnitudes and colours that 
place them within the RGB locus shown in Fig. 3. Because 
our sample of 44 RGB candidates constitutes only a small 
fraction (« 13%) of the full sample of stars in this region 
of the CMD and lying within ~ 250"of the cluster core, we 
conclude that we would expect, at most, 1 interloper in our 
final sample of 40 cluster members. 

We now proceed by calculating the average radial veloc- 
ity and velocity dispersion of NGC 2419 using the method of 
Pryor & Meylan (1993). We first assume that each velocity, 
v iy is drawn from a normal distribution, 

m = , 1 ™p (- prrt^ ) > ( 3 ) 

where v c and a c are the cluster's average radial velocity and 
intrinsic velocity dispersion, and a eA is the individual veloc- 
ity error for each star. Calculating the likelihood function, 
L, for all stars and taking the partial derivatives of its log- 
arithm I = logL with respect to v c and a c leads to the 
following set of equations: 

N N 

E ? 2 r 2 x ~^E rrj-in = ( 4 ) 
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The above equations can be solved analytically to obtain a 
first guess value for v c and a c if zero measurement errors 
a e< i for all stars are assumed. The full solution can then 
be obtained iteratively, starting with the solution for zero 
measurement errors as first estimate. The errors of v c and 
<T C are obtained from the information matrix, /, by 



(6) 
(7) 



where the components Iij of the information matrix are cal- 
culated as given in Pryor & Meylan (1993). 

Using the 40 probable members from Table 1, we obtain 
a mean cluster velocity of v c = —20.63 ± 0.74 km/sec and 
an intrinsic velocity dispersion of a c = 4.61 ± 0.53 km/sec. 
This mean velocity agrees very well with the one obtained 
by Olszewski, Pryor & Schommer (1993). Our velocity dis- 
persion is, however, significantly higher than the value of 
2.7 ± 0.8 km/sec determined by Olszewski, Pryor & Schom- 
mer (1993) and agrees with their value only at the 2a level. 
The explanation for this difference appears to lie in the ra- 
dial distribution of the stars in the two samples: Fig 5 reveals 
the Olszewski, Pryor & Schommer (1993) sample to contain 
only three stars in the central r ~ 1', where we observe a 
significant rise in the cluster velocity dispersion based on a 
larger sample of 17 stars in this region. In the outer parts of 
the cluster, our measured velocity dispersion (see below) is 
in very good agreement with the previous estimate. 

The lower panel of Fig. 5 shows the distribution of radial 
velocities of the 40 candidate members in the velocity range 
-36 to -4 km/sec (open histogram). The dotted curve shows 
a gaussian fit to this distribution using the above values 
for v c and a c . As can be seen, the observed radial velocity 
distribution is well approximated by this Gaussian except, as 
noted above, that the number of outliers in the wings of the 
distribution is somewhat higher than predicted: i.e., stars 
S17 and S38 are more than 2.2<r away from the mean while 
star S22 is almost 3cr away from the mean. For a Gaussian 
distribution and a sample of 40 stars, one would expect to 
find a star that is more than 3<r away from the mean in only 
10% of all cases. Also, the chance to have three or more stars 
more than 2.2<r away from the mean is only 10%, while the 
chance to have two such stars is twice has high. This makes 
it possible that at least one star, most likely star S22, is 
either a cluster nonmember or a radial velocity variable. 

Since star S22 is located in the cluster center and has a 
photometry that places it squarely on the RGB (see the right 
panel of Fig. 3), it could be a binary system, in which case 
it should be excluded from the analysis. The absolute radial 
velocity difference between star S22 and the cluster mean 
(~15 km/sec) is also so large that star S22 would be nearly 
unbound if the radial velocity difference is due to a different 
orbital velocity. Star S22 is also one of the five stars in com- 
mon with Olszewski, Pryor & Schommer (1993), and has the 
most discrepant velocity in both samples. Moreover, two in- 
dependent measurements different by 5.1 km/sec, or three 
times the quadrature sum of the individual uncertainties. 
These facts suggest that star S22 may be a cluster binary 
or, perhaps more likely, an RGB star that shows a velocity 
"jitter" that has previously been found for some globular 
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Figure 5. Upper panel: Heliocentric radial velocity for NGC 2419 
probable members plotted as a function of distance from the clus- 
ter center. The dashed line shows our best estimate for the cluster 
systemic velocity, (v r ) ~ —20.3 km/sec. There is a clear rise in 
velocity dispersion toward the cluster center. Note the discrepant 
measurements for star S22, which is the largest outlier in both our 
sample and that of Olszewski, Pryor & Schommer (1993). Lower 
panel: Radial velocity distribution of the measured stars in NGC 
2419. Dotted and dashed lines show Gaussian fits to this distribu- 
tion obtained by either using all stars from Table 1 (dotted line), 
or by excluding star S22 — which has the most discrepant ve- 
locity — from the sample (dashed line). The filled red histogram 
shows the expected distribution, in this velocity interval, for fore- 
ground stars within a 500" X 500" region along this line of sight 
and falling along the RGB in the cluster CMD (sec the right panel 
of Fig. 3). The solid line shows the distribution of stars from our 
model cluster fitted to the data without star S22. It follows very 
closely the corresponding Gaussian distribution (dashed line). 



cluster stars on the upper RGB (e.g., Gunn & Griffin 1979, 
Mayor et al. 1984). If we omit star S22, we obtain a mean 
cluster velocity and an intrinsic velocity dispersion of: 



-20.28 ±0.68 
4.14 ±0.48 



km/sec 
km/sec. 



(8) 



These values are within the errorbars of the values that we 
obtain when using all stars and we will use them throughout 
the rest of this paper. 

We note that the three stars with the most discrepant 
radial velocities are also among the four most central stars. 
While this could simply be a statistical effect, it could also 
point to either a high binary fraction in the core or the 
broadening of the velocity distribution due to additional un- 
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Figure 6. Radial velocities as a function of position angle (PA) 
measured from north (PA=0°) towards east. The solid curve 
shows our best-fitting rotation curve for NGC 2419. It has an 
amplitude of 3.26 ± 0.85 km/sec, implying significant rotation 
of the cluster. The direction of rotation, PA=40.9° ± 17.8° is, 
within the errobars 90° away from the semi-major axis found in 
the surface density profile by Bellazzini (2007), PA=+105° ± 28° , 
supporting the fact that NGC 2419 is rotating. 



seen mass, possible if the cluster is mass segregated or con- 
tains an intermediate-mass black hole (Baumgardt, Makino 
& Hut 2005). Multi-epoch observations of the stars in our 
sample and radial velocities of additional core stars would 
be required to test these interesting possibilities. 



3.3 Cluster Rotation 

Our sample of 40 member stars also allows us to check for 
a possible rotation of NGC 2419. Fig. 6 shows the radial 
velocities of stars as a function of position angle (PA). Here 
the position angle is measured from north (PA=0°) towards 
east (PA=+90°). There is a clear dependence of the average 
velocity from the position angle visible in the data. This is 
confirmed by fitting a sinusoidal curve to the radial velocity 
data. Allowing for a variable rotation angle, we obtain a ro- 
tation amplitude of 3.26 ± 0.85 km/sec around an axis with 
position angle of PA = 40.9 ± 17.8 degree. Our data there- 
fore implies cluster rotation at more than the 3<r level. The 
root mean square scatter around the rotation curve is 4.00 
km/sec. NGC 2419 is therefore partly rotationally supported 
and partly supported by random stellar motions. Bellazzini 
(2007) found that NGC 2419 is slightly elliptical with aver- 
age ellipticity of e = 0.19 ± 0.15 and position angle of PA 
= +105° ± 28°. Both the relatively small amount of ellip- 
ticity and the average angle of ellipticity, which is within 
the errorbars 90° away from the rotation axis that we find, 
supports our finding that NGC 2419 is rotating. 



3.4 Mass modeling and mass-to-light ratio 

The global mass-to-light ratio is calculated from the veloc- 
ity dispersion and the total cluster luminosity under the 
assumption that mass follows light. Unless there is primor- 
dial mass segregation, this is probably a valid assumption 
since the relaxation time of NGC 2419 is larger than a Hub- 
ble time: i.e., the central and half-mass relaxation times are 
T r o « 10.5 Gyr and T r h ~ 19 Gyr, respectively (Djorgovski 
1993). In addition, Dalessandro et al. (2008) found no ev- 
idence for mass segregation among the blue straggler stars 
in NGC 2419, which supports the assumption that mass fol- 
lows light. We also neglect the cluster rotation and assume 
that the cluster has an isotropic velocity dispersion. In such a 
case, we can use the method described in Hilker et al. (2007) 
to calculate the mass-to-light ratio. Their method first de- 
projects the surface density profile (here a cored Sersic pro- 
file) to obtain the three-dimensional spatial density distri- 
bution, p(r). The spatial density is then used to calculate 
the potential, <E>(r), of NGC 2419 and from p and the dis- 
tribution function, f(E), is calculated by assuming spheri- 
cal symmetry and using eq. 4- 140a of Binney & Tremaine 
(1987). This distribution function is then used to create an 
A-body representation (i.e., particle positions and veloci- 
ties) of NGC 2419 with a given first guess of M' for the 
cluster mass. We then draw a number of N' stars located 
at the same projected radius of each observed star from this 
A-body model and measure the projected velocity disper- 
sion, <TMod, of the sample stars. The total mass of NGC 2419 
can then be calculated by comparing this velocity dispersion 
with the observed velocity dispersion according to: 

M c = M' ■ a 2 Mod /a 2 c . (9) 

This method has the advantage that it can work for any 
given observed surface density profile and for any given ra- 
dial distribution of stars. This sets it apart from parameter- 
ized formulas which work only for certain density distribu- 
tions and usually need the core or global velocity dispersion 
to calculate the cluster mass. If the stars with radial ve- 
locity measurements are neither concentrated in the cluster 
core, nor spread out over the cluster in the same way as the 
cluster stars, using such formulae can lead to a bias in the 
derived cluster mass. 

We calculate an A-body model with A = 5 • 10 stars 
in total and extract for each observed star N' — 200 stars 
from this model. Calculating the velocity dispersion for the 
theoretical model and comparing it with NGC 2419, we find 
a total mass of Ma = (9.02 ± 2.22) • 10 5 M for NGC 2419 
and a mass-to-light ratio of M/Ly — 2.05 ± 0.50 in solar 
units, using all stars except star S22. If we were to include 
star S22, the total mass would rise to M c = (1.09±0.26)-10 6 
Mq and the mass-to-light ratio would be M/Ly = 2.48 ± 
0.60 in solar units. These values are significantly higher than 
the M/L ratio derived by Olszewski et al. (1993), M / L v ~ 
0.7 ± 0.4 Mq/ Lv,q- Since our method gives a very similar 
M/L than what they found if we use their stars as input, 
the main reason for the discrepant M/L values seems to be 
the different input stars. 

We finally tested our modeling method by randomly 
drawing 40 stars from the A^-body model and then measur- 
ing their velocity dispersion and fitting the model against 
them in the same way as was done for the NGC 2419 data. 
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Figure 7. Observed velocity dispersion as a function of radius. 
The solid line shows our prediction based on the best-fitting cored 
Sersic model. Dotted and dashed lines show the predicted velocity 
dispersion if we add NFW halos with a scale radius of Rg = 500 
pc and masses of M DM = 4 ■ 10 6 M and M DM = 10 7 M Q 
inside Rg to this model. Models with additional dark matter halos 
significantly overpredict the velocity dispersion in the outer parts, 
showing that NGC 2419 does not possess a dSph-like DM halo. 

Doing this 100 times, each time with a different random re- 
alisation of stars, we obtained an average ratio of derived 
mass-to-light ratio to true mass-to-light ratio of 1.00 ±0.01, 
i.e. our modeling method gives an unbiased estimate of the 
true mass-to-light ratio. 

In order to compare the observed mass-to-light ratio 
with the predictions of population synthesis models, we need 
to know the cluster age and metallicity. Based on an analysis 
of a high-resolution Echelle spectrum for a single RGB star 
(10), Shetrone, Cote & Sargent (2001) measured a metallic- 
ity of [Fe/H]=-2.32 for NGC 2419, which is one of the lowest 
abundances among all Galactic globular clusters. Salaris & 
Weiss (2002) have used horizontal branch and turn off star 
brightnesses to determine ages of Galactic globular clusters 
and derived an absolute age of 12.3 Gyr or 12.8 Gyr for 
NGC 2419, depending on the adopted metallicity scale. For 
this metallicity and age, the SSP models of Maraston (2005) 
predict a mass-to-light ratio of 1.95 for a Kroupa (2001) 
IMF, while the Bruzual & Chariot (2003) models predict 
M/L v = 1.85 for a Chabrier (2003) IMF. These values agree 
remarkably well with the derived mass-to-light ratio. The 
M/Lv ratio of NGC 2419 is therefore entirely compatible 
with the assumption that the cluster formed with a stan- 



dard stellar mass function. This sets NGC 2419 apart from 
some UCDs, which have only slightly larger sizes and masses, 
and for which standard mass functions seem to underpredict 
their mass-to-light ratios (Hasegan et al. 2005; Mieske et al. 
2008). 

The "normal" mass-to-light ratio of NGC 2419 also ar- 
gues against the presence of a significant amount of dark 
matter in the cluster. In order to test if a possible DM halo 
could exist at least in the outer cluster parts, we have di- 
vided the sample into three equally large groups according 
to their radial distance and calculated velocity dispersions 
separately for each group. Fig. 7 compares velocity disper- 
sions at different radii with the predicted velocity dispersion 
profile from our iV-body model. It can be seen that the ve- 
locity dispersion of the cored Sersic model is in reasonable 
agreement with the data. If at all, additional mass is required 
only in the inner parts of the cluster, where the predicted 
velocity dispersion is smaller than the observed one by about 
la, not in the outer parts as might be expected for a dark 
matter halo. 

We also added different DM halos to our iV-body model 
in order to see the effect of these halos on the velocity dis- 
persion. The DM halos follow NFW profiles with scale radii 
Rs ~ 500 pc, which agrees with measured scale radii of 
Galactic dSph galaxies which are generally larger than a few 
hundred pc (Wu 2007), and had various ratios of dark halo 
mass inside Rs to the stellar mass of NGC 2419. We then 
calculated a new velocity dispersion profile for NGC 2419 
and fitted it to the data in the same way as described above. 
From our fit procedure, we derive a best-fitting (stellar) M/L 
ratio and the total stellar and dark mass of NGC 2419. As 
can be seen in Fig. 7, significant amounts of dark matter 
are basically ruled out by our observations since they over- 
predict the velocity dispersion in the outer cluster parts, at 
least for models with isotropic stellar velocity dispersions. A 
DM halo with total mass of Mdm = 1-65 ■ Mstars overpre- 
dicts the velocity dispersion by more than 2a at a radius of 
R=140". This model also has a best-fitting stellar M/L ratio 
of only M/L = 1.40±0.34, which is 2a below the predictions 
of the stellar evolution models. Our observations therefore 
rule out a DM halo of more than M « 1 • 10 7 M Q inside 500 
pc around NGC 2419, corresponding to a density of 0.02 
Mq/pc 3 . This is significantly lower than the dark matter 
content of Galactic dSph galaxies inside a similar radius as 
derived by Walker et al. (2007) and Strigari et al. (2008). It 
is also a factor 5 lower than central dark matter densities of 
dwarf galaxies as determined by Gilmore et al. (2007). Our 
observations therefore argue strongly against the formation 
of NGC 2419 within a sizeable dark matter halo. More ra- 
dial velocities would help to strengthen our conclusions and 
test the validity of several assumptions that we made for 
our data analysis, in particular that the stellar orbits are 
isotropic and that the cluster is not mass segregated. 

Finally, we note that the low observed velocity disper- 
sion in the outer parts of NGC 2419 might also be a prob- 
lem for MOND. According to MOND, the orbits of stars 
show a difference with respect to standard Newtonian be- 
havior once their acceleration falls below a critical acceler- 
ation a ~ 1.2 • 10~ 8 cm/sec 2 (Sanders & McGaugh 2002). 
While the acceleration of stars near the half-mass radius of 
NGC 2419 is still higher than ao, it falls significantly below 
do in the outer parts of NGC 2419. At the radius of our out- 
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crmost datapoint in Fig. 7, the internal acceleration is only 
one third of ao, while the external acceleration is only one 
tenth of ao, so one would expect that the observed veloci- 
ties start to deviate significantly from the prediction of our 
iV-body model which was created based on Newtonian dy- 
namics. This is not seen in our data. In fact, as shown by Mil- 
grom (1994), if one neglects the external field of the Milky 
Way, MOND predicts that the line-of-sight velocity disper- 
sion should level off at a value a m in = 0.471-^/ GMc/ao at 
large distances. If we use the above calculated cluster mass, 
we find a min = 2.6 km/sec for NGC 2419. Our outermost 
datapoint is already probing this value, so acquiring addi- 
tional data at even larger radii might prove a powerful way 
to test the validity of MOND. 



4 CONCLUSIONS 

We have measured precise radial velocities for w 40 stars in 
the outer halo globular cluster NGC 2419 using the Keck 
telescope High Resolution Echelle Spectrometer. We find a 
mean cluster velocity of <v r >— —20.3 ± 0.7 km/sec and 
an internal velocity dispersion of a = 4.14 ± 0.48 km/sec. 
Our observations also reveal a slight cluster rotation with 
amplitude 3.26 ± 0.85 km/sec around a position angle of 
40.9 ± 17.8 degrees. From a comparison of the measured 
velocity dispersion to the velocity dispersion of a spherically 
symmetric model for NGC 2419, we find a total mass of 
M c = (9.01 ± 2.22) • 10 5 M and a mass-to-light ratio of 
M/L v = 2.05 ± 0.50 in solar units. This value is entirely 
compatible with the one expected for a stellar system at the 
metallicity of NGC 2419 following a standard mass function 
like Kroupa (2001) or Chabrier (2003). 

NGC 2419 therefore does not show any dynamical ev- 
idence of a significant depletion of low-mass stars, which is 
consistent with expectations given the large relaxation and 
dissolution time for this cluster. Similar M/L ratios which 
are in agreement with standard mass functions have been 
found for other galactic globular clusters with large relax- 
ation and dissolution times like Omega Cen (M/Ly = 2.5 
van de Ven et al. (2006)). This sets Milky Way globular clus- 
ters apart from UCDs for which standard mass functions on 
average seem to underpredict mass-to-light ratios (Mieske 
et al. 2008). If real, this points to a significant change in the 
star formation process occurring at a characteristic mass of 
M c « 2 • 10 6 M Q (e.g., Haseegan et al. 2005; Mieske et al. 
2008, Dabringhausen et al. 2009). Such a transition might 
arise, for example, if more massive clusters become optically 
thick to far infrared radiation and are born with top heavy 
initial mass functions (Murray 2008). 

We do not find any evidence for the presence of sub- 
stantial amounts of dark matter in NGC 2419. Since both 
the depletion of dark matter from the cluster center due to 
two-body relaxation and dynamical friction (Baumgardt & 
Mieske 2008), as well as tidal stripping of dark matter from 
the cluster halo are unlikely, this indicates that NGC 2419 
did not contain substantial amounts of dark matter at the 
time of formation. NGC 2419 is, however, one of the most 
likely candidates for a globular cluster to have formed with 
an associated dark matter halo given to its quite low metal- 
licity. Our non-detection of dark matter in this cluster there- 
fore supports the hypothesis that globular clusters, as a rule, 



did not form at the centers of dark matter halos. Instead, 
an origin driven by gas dynamical processes during merg- 
ers between galaxies (Bournard, Due & Emsellem 2008) or 
proto-galactic fragments seems to be the more likely expla- 
nation for the formation of even the lowest metallicity glob- 
ular clusters. 

We also find a slightly larger than expected velocity 
dispersion in the central regions of the cluster, which might 
point to a significant binary fraction or additional unseen 
matter in the core. At the same time, the rather low veloc- 
ity dispersion in the outer regions of the cluster could pose 
a problem for alternative gravitational theories like MOND. 
These conclusions rely on several assumptions that we had 
to make for our modeling, in particular that stellar orbits are 
isotropic throughout the cluster and that the stellar mass- 
to-light ratio is constant with radius. Acquiring additional 
radial velocities would help test the validity of these assump- 
tions and strengthen our conclusions. 
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